# libraries
librarian::shelf(
here, mapview, sf,
quiet = T)
# paths
tri_shp <- "/Users/bbest/My Drive/projects/bbnj/data/polynesian-triangle_shapefile/polynesian_triangle_clipped.shp"
# read
tri <- read_sf(tri_shp)
tri
Simple feature collection with 1 feature and 1 field
Geometry type: MULTIPOLYGON
Dimension: XY
Bounding box: xmin: 1828474 ymin: -5958241 xmax: 11217990 ymax: 3177202
Projected CRS: WGS 84 / PDC Mercator
# A tibble: 1 × 2
FID geometry
<dbl> <MULTIPOLYGON [m]>
1 0 (((4295005 -5288629, 4294399 -5286727, 4293674 -5284491, 4292939 -52822…